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Reduction-Based Robustness Analysis of Linear 
Predictor Feedback for Distributed Input Delays 

Anton Ponomarev 

Abstract —Lyapunov-Krasovskii approach is applied to parameter- and 
delay-robustness analysis of the feedback suggested by Manitius and 
Olbrot for a linear time-invariant system with distributed input delay. 
A functional is designed based on Artstein’s system reduction technique. 
It depends on the norms of the reduction-transformed plant state and 
original actuator state. The functional is used to prove that the feedback 
is stabilizing when there is a slight mismatch in the system matrices and 
delay values between the plant and controller. 

Index Terms —Delay systems, predictive control for linear systems, 
robust control, Lyapunov methods. 

I. Notation 

We write M > 0 or M > 0 to state that a symmetric real matrix 
M is positive definite or positive semidefinite, respectively. Also in 
this case A m i n (M) and A ma x(M) represent the minimal and maximal 
eigenvalues of M. Vector norms being used are ||a;|| = V x T x and 
||s;||m = Vx T Mx, where M > 0. Euclidean matrix norm is ||M||. 

The symbol PC(T,X) stands for the space of piecewise contin¬ 
uous functions mapping T C R into a Euclidean space X. The L 2 
norm of <p G PC([—h, 0), f? r ) is ||<p||, i.e., 

IMI 2 = f M0)\\ 2 de. (i) 

J-h 

Given u G PC(R, R r ), let ut be a function defined as Ut{6) = 
u(t + 9) for all 9 G [—h, 0). The constant h is specified below. 

II. Introduction 

A. The problem 

Consider the time-invariant system 

n ,o 

x(t) = Ax(t) + ^ Biu(t — hi) + / B in ,(9)u(t + 9) d9, (2) 

i= 1 ■'-him 

where x G R n , u G R r , h\ > 0, /lint > 0, and £>i n , G 
PC([-h int , 0],f? nxr ). For brevity, we will use Stieltjes integral 
notation and write the system under consideration as 

x(t) = Ax(t) + f d/3(9)u(t + 9), (3) 

J-h 

where h > max{/u, / 12 ,..., /uv, /lint}, 

13(9) = J2Bi X {9 + hi)+ / Bi„, (t) dr, (4) 

i= 1 “ hint 

and x is the Heaviside step function. 

The following control law was proposed for 0 in (TJ: 

/ 0 rO 

J e A(T ~ e) d(3(r)u(t + 9) d9, (5) 

where F is a constant matrix. The feedback 0 is called a predictor 
feedback because it employs the plant’s model (i.e., the matrices A 
and /3) to, in a sense, predict the future state of the plant. 

Our goal is to investigate robustness of the feedback 0. In terms 
of 0. we are interested in: 

• parametric robustness (small uncertainty in A, Bi , and B i„,); 

• delay-robustness (small uncertainty in hi and /i; nt ). 
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B. Previous results overview 

A range of methods is known to be suitable for analysis of linear 
systems of the form 0. 0. Let us separate them into those using 
Lyapunov-Krasovskii functional analysis and those doing otherwise. 

Most of the progress with non-Lyapunov techniques has been 
achieved in the area of systems with one discrete delay, e.g., delay- 
robustness of a predictive controller 0, 0, robustness with respect 
to a finite-sum implementation 0, 0, robustness of an adaptive 
controller in presence of a disturbance (6), and delay-robustness of 
a linear time-varying predictor feedback 0- Furthermore, it has 
been shown in 0 that robustness with respect to a finite-sum 
implementation may be ensured by including a low-pass element in 
the control loop. 

Lyapunov-Krasovskii analysis of systems with one discrete delay 
was shown to succeed in proving delay-robustness of the predictor 
feedback (9) and robustness with respect to uncertain parameters 
ED- Adaptive controllers were designed in HD, DU- Recently, a 
predictor feedback for retarded GUI and neutral m systems with 
state delays and an input delay was proposed, the closed loop’s 
exponential stability being proven with a functional as well. 

Lyapunov-Krasovskii analysis has been performed for distributed 
delays too but less extensively. The results are closed-loop ex¬ 
ponential stability (Hi and stability with respect to an additive 
disturbance GU. This paper expands the list with parameter- and 
delay-robustness of the feedback. 

C. Summary of the note 

In Section III. the loop 0. 0 is turned into 0, 0 using the 
transformation 0 borrowed from GU- The Lyapunov-Krasovskii 
functional G3 is then constructed. It includes the norms of the 
transformed plant state y(t) and original actuator state ut . Lemmas 
00 prove that this functional is quadratically bounded. 

In Section IV, we show that a mismatch in A and j3 between 
0 and 0 introduces a distributed delay into the controller part of 
otherwise delay-free transformed system 0. 0 which now becomes 
0. 03i. The delay does not significantly affect the system behavior 
if the mismatch is negligible. It leads to the main result : closed-loop 
stability is robust (Theorem 0. In order to facilitate a comparison 
of our approach with the preceding ones, we provide Corollary 0 
together with its concise proof which is the special case of Theorem 
0for systems with one discrete delay. 

III. Lyapunov functional 

Following Artstein GU, let us introduce a new variable 

y(t) = x(t) + f Q(9)u(t + 9) d0, (6) 

J-h 

where 

r e 

Q(9)= / e A(T ~ e) dp(r). (7) 

J-h 

The closed loop 0,0 in the new variables takes the form 

y(t) = Ay(t) + Q(0)u(t), (8) 

u(t) = Fy(t). (9) 

It has been demonstrated by direct calculation in 0 that the 
eigenvalues of the closed loop 0, 0 coincide with the eigenvalues 
of the matrix A + Q(0)F. Suppose the matrix is Hurwitz, so that the 
nominal closed loop is exponentially stable. 

.00 © 2015 TF. F.F. 
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To apply Lyapunov approach, we aim at finding a functional 
v(x,p) defined for all x £ R n and ip £ PC(J—h, 0), f? r ) that 
admits upper and lower bounds proportional to ||ai|| 2 + ||<p|| 2 . We 
first construct a functional of this kind for the transformed loop GO, 
© and then come back to the original variables. 

Let us choose arbitrary matrices W' > 0 and W" > 0. Suppose 
then that V > 0 is the solution of 


{A + Q(0)F) t V + V{A + Q(O)F) = ~{W' + 2 F t W"F). (10) 
For the loop GQ, we propose the functional 


v(y,p) = \\y\\v + 


e° e \\pm\w" d9 


defined for all y £ R n and tp £ PC([—h, 0), R r 'j , where 

A min (W') 


A max (V) ’ 


( 11 ) 


( 12 ) 


In the original variables it is 
r o 


v(x,p) = 


x+ Q(e)p(e)de 
J-h 


+ 


e M0)ll W"d0. (13) 


The following two lemmas prove that the functional ( 1131 has the 
required properties (upper and lower bounds). 

Lemma 1: 

v(x,p) < M(\\x\\ 2 + \\<p\\ 2 ), (14) 

where 


M = max {2A max (V) ||e 


Ah 11 2 


A max (LL") + 2/iA max (l/) max ||Q(0)|| 2 }. (15) 

9e[0,h] J 

Proof: One obtains this from GEU using the Young’s inequality, 
the triangle inequality, and the Cauchy-Schwarz inequality. ■ 

Lemma 2: 


v(x, tp) > m-u 


v(x, i p) > m a; ||a;|| 2 , 


(16) 


where 


rriu = e ^Amin (W"), 

(17) 

m x — A m i n ^ ^ V + m u G) ^, 

(18) 

G= f Q{0)Q T (0) d0. 

Jo 

(19) 


Before we proceed to the proof, take notice that the lower bound 
of the functional in terms of the full norm follows from Lemma [2] 

v(x,ip)> --- ( 11*11 +MI )• ( 20 ) 

We leave the lemma as it is, though, lest the estimations be more 
conservative than necessary. 

Proof: The initial step is to estimate 

v(x,ip) > \\y\\v + m u ||v5|| 2 , (21) 

where y is linked to x and tp via ©: 

y = x+ f Q(0)ip(0)d9. (22) 

J-h 

The first inequality is obtained by dropping ||y||y in ( 1211 . 

The idea used to establish the second inequality is to consider 
an optimal control problem: minimize the right hand side of Oil 
with respect to tp. Any function ip £ PC([—h,0),R r ) allows the 
decomposition (orthogonal projection on the rows of Q) 

<p(9)=Q T (6)c + f J (6), (23) 


where c is a constant vector and ip £ PC(\—h, 0), _R r ) satisfies 


Q(9)ip(9)d9 = 0. 


(24) 


This representation lets one write 


y = x + Gc, (25) 

\\ p \\ 2 = c T Gc + \\ il )\\ 2 > c t Gc (26) 

=> v(x, ip) > ||Gc ||v + m u c T Gc + 2c T GVx + ||*||v. (27) 


The minimum of the quadratic estimation is at 

c = — (G + m„T _1 ) 1 x. (28) 

We substitute this c into ( 1271) to find that 

v(x,ip) > * T (V _1 + m a 1 G) X x, (29) 

which leads to the desired inequality. ■ 


IV. Robustness analysis 


A. General result 

The controller’s robustness with respect to a mismatch in the 
prediction model is analyzed here. In this scenario, the exact plant 
model 0 or 0 is unknown but its estimation 

N ,0 

x(t) — Ax(t ) + Bju(t — hi) + / B m i{9)u(t + 9) d9 (30) 

i= 1 J -hint 


or, equivalently, 


x(t) = Ax(t) + / dj3(9)u(t + 9) (31) 

J-h 

is available. Here hi > 0, h. mt > 0, B i„t £ PC(j— h m t, 0], R nxr ^j , 

h > max{/ii, fi 2 , • ■ ■, h.N, hmt}, 

N ^ ^ r m a x{S,-ii int ) 

0(6) = y' j Bix{9 + hj) + / Bmt(r)dT, (32) 

i=l d-h m , 


and x is the Heaviside step function. Observe that h is the same 
in the nominal system 0 and its approximation 0}: for that, it is 
sufficient to take 


h ^ max{/ii, , • ■ •, Itiv, hmt , h\ , /t 2 , ■ ■ ■, Itiv, Itint (■ 

The controller designed from ( 130b or ( 131b would be 

u(t) = F (x(t) + [ Q(9)u(t + 9) dd\ , 


where 


Q(0) = 


AG-6) 


dP(r). 


(33) 

(34) 

(35) 


Regarding the choice of F, we demand that A + Q(0)F be Hurwitz. 
However, matrices A and Q(0) are not known exactly due to para¬ 
metric uncertainties. Nevertheless, suppose that one may establish 
some boundaries on A and Q(0) and choose a value of F which 
renders A + Q(0)F Hurwitz for all possible values of A and Q(0). 
This problem is not in the focus of the paper, so we take such F as 
a given and assume hereafter that A + Q(0)F is Hurwitz indeed. 

After the transformation 0, the plant is still 0 and controller 
03 is written as 


u(t) = F ( y(t) + J AQ(9)u(t + 9) d9 


(36) 
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where A Q(9) = Q(0) — Q{9). One can clearly see how control 
delay reappears in the transformed loop 0. 1361 due to imperfect 
modeling. 

Let v(t) be the value that the functional 113b takes on a specific 
solution of the closed loop 10. 134b . 

Lemma 3: Along the solutions of the closed loop 0. 134b the 
functional G3 satisfies 

v(t ) < —ov(t), (37) 


where 


a = a — fci||AQ|| — fe||AQ|| , 
\\AQf=[ \\AQ(d)\\ 2 d6, 

J-h 


ki = 


l|VQ(0)F|| 


k 2 = 


2Amax(W")||F|| 2 


min{Amin(f / ),m u } ’ 
a is 112b . and m u comes from Lemma [2] 


m u 


Proof: Differentiating v(t ), we use 110b to get 
v{t) < 2y T (t)V(Ay{t) + Q(0)u(t)) 

r° 

— <j I e \\u(t + 9) ||vu /7 d9 + ^(f)live" 


(38) 

(39) 

(40) 


(41) 


< -ov(t) + 2 y 1 (t)VQ(0)F / A Q(9)u(t + 9) d9 

J-h 

+ 2||w(i) — Fy{t)\\w" ■ 


Further estimations include 
rO 


and 


2y 1 (t)VQ(0)F AQ(9)u(t + 9) d9 

J-h 

< 2||HQ(0)F||||AQ||||t / (t)|||M| 

< ||yQ(0)F||||AQ||(||t/(f)|| 2 + || Mt || 2 ) 

< fci||AQ|| v(t) 

2\\u(t) - Fy(t)\\ 2 w „ <2A ma x(fL")l| J F|| 2 (l|AQ|||| Wf ||) 5 
< k 2 \\AQ\\ 2 v{t). 


(42) 


(43) 

(44) 

(45) 

(46) 

(47) 


In the end, we arrive at the estimation v(t) < —av(t). ■ 

The following is our main result. 

Theorem 1: If A + Q(0)F is Hurwitz and 130b approximates 0 
closely enough that (see the remark after the proof) 

n2 | s/kf + 4:k2a — ki \ 


iiAoir < 


2 k 2 


(48) 


Remark: Let us explain why ||AQ|| becomes small when approx¬ 
imation of 0 with 130b is almost perfect. If all parameters of 130b 
closely match those of 0. then (3 and (3 by definition are uniformly 
close on [— h, 0] except for a finite number of small intervals like 
[-hi,-hi), [-hi,-hi), [-hi nt ,-h ml ) or [-h mt , -/lint). It follows 
then that AQ is uniformly small on [— h, 0] except the small intervals 
where it is bounded (and the bound is independent of the intervals’ 
lengths). Smallness of ||AQ|| ensues. 

See the proof of Corollary |T| for a quantitative example of how 
|| AQ|| depends on a mismatch in the value of a discrete delay. 


B. Special case: one discrete delay 

To highlight the conditions imposed by Theorem |T| on discrete 
delay values, we supplement it with an application to the widely 
studied single discrete delay case. 

Corollary 1: The closed loop 


x(t) = Ax(t) + Bu(t — 5), 
u(t) = F ^ 'x(t) + J ~ e~ A ^ +e) Bu(t + 9) d9 


(56) 

(57) 


where S > 0, 5 > 0, is exponentially stable if there exist V > 0, 
W' > 0, W" > 0 such that 

(.A + e~ A5 BF) T V + V(A + e~ AS BF ) =-{w' + 2 F t W”F) 

(58) 

and 


|5 ” 5|< p| ln(l + 


s/l + k 3 e 2 \\ A \\ h - 1 
e 2||A ||h 


where 


kz = 


2||A|| / \Jk\ 4- 4fc 2 cr - ki 


||B|| 


ki = — 


2 k 2 

\\Ve~ AS BF\\ 


min (A min (y), e—^A min (TL")} ’ 
2A ma x(fL")||F|| 2 


k2 e— h \ min (W") ’ 

A min (W / ') t r , ?1 

erm— - h = max{6,5}. 

^max \ V J 

If A = 0, then the delay mismatch condition is 


|S-5| < 


\Jk\ + Ak 2 cr — ki 

2||B||fc 2 


(59) 

(60) 

(61) 

(62) 

(63) 

(64) 


then the closed loop 0, (134b is exponentially stable: 



l*(f)H 2 < e (||*(0) || 2 + ||mo|| 2 ), 

11 Lx 

(49) 


l|Mi|| 2 < — e -<Tt (ll*(0)|| 2 + llnof). 

TYlu 

(50) 

Here || AQ||, ki, k 2 , M, m x , m u , and a are defined in Lemmas00 

and we mean 

that (130b approximates 0 in the sense that 


A 

~ A, 

(51) 

Bi 

w Bi for all* = 1, 2, ..., N, 

(52) 

Bmt 

S3 Bint uniformly on [ — min{/imt, /lint}, 0], 

(53) 

hi 

w hi for all * = 1, 2,..., N, 

(54) 

hint 

?3 /li n t. 

(55) 

Proof: It follows from Lemmas Q}0 

■ 


Proof: Given a solution ( x(t),Ut j of the closed loop 156b . 157b . 
we study the behavior of 


r° 

)-h 


v (t) = l|y(*)llv + / e CT8 ||w(f + 0)||^" d0, 


where 


y(t) = x(t) + J Q(0)u(t + 0) dO, 
v \o, if 6* < —5. 

The control law ( 157b is equivalently represented as 

u(t) = F (y(t) + [ AQ(0)u(t + 0) dO 


(65) 

( 66 ) 

(67) 

( 68 ) 
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where A Q(0) = Q(0) — Q(0) and 


References 


Q{6) = 


e~ A(S+e) B , if (9 e [-5,0], 

0 , if e < -8. 


(69) 


Differentiating v(t) and performing some estimations explained in 
the proof of Lemma [3] we get 

v(t) < (-a + fci||AQ|| +fc 2 ||AQ|| 2 )v(f), (70) 


where 


If 


l|AQ|| 2 = / \\AQ(e)f dd. 


-ff + fci||AQ|| + fe,||AQ|| 2 < 0, 


(71) 


(72) 


then exponential stability of d56L d57b is guaranteed by Lemmas[T|{3] 
Let us estimate ||AQ||. If h = 8 > 8, 


A Q(9) = 




if 9e [—h, —8), 


(j - e A ^~ s ^ e - A< ~ S+e) B, if 9 e [-5, 0]. 


(73) 


||AQII > s small because A Q(9) is bounded on the first interval 
which has small length, and small on the second one which length 
is bounded: 

r-s r o 


l|AQ|| = / \\AQ(9)\\ 2 d9 + / ||A Q(9)\\ 2 d9 


< 


\\ B \\ 

2pi 


+ 


2||A|| 


l«-«l 


2 \\ A\\h 


-1 


■>(< 


|i-«l 


-1 


If A = 0, then by continuity d75b turns into 
IIAQII 2 < ||73|| 2 |5 — <5|. 

If h = 8> 8, 

e~ A(h+e) B, if 9 € [~h, -8), 

(I - e- A| ^- 5|S ) e~ A(S+e) B, if 9 e [-5,0], 


A Q(9) = 


(74) 


(75) 


(76) 


(77) 


but the same estimation GD or {76} holds true. The corollary’s 
premise then resolves in {72}, thus proving exponential stability. ■ 
It should be mentioned that Corollary |T| is not the best result 
available for the single discrete delay case as our estimations are 
quite conservative. The problem has been widely studied with more 
accurate results achieved, e.g., in £1, 0. 
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V. Conclusions 

A Lyapunov-Krasovskii functional has been devised for a linear 
time-invariant system with distributed input delay closed by a predic¬ 
tor feedback. The functional is based on Artstein’s system reduction 
technique. 

The framework is shown to be convenient for robustness analysis 
of the closed loop: it alleviates the proof of closed-loop exponential 
stability when controller’s predictive model (matrices and delay 
values) differs slightly from the actual plant. 

Future research may include a generalization of the approach for 
nonlinear systems and systems with both state and control delays. 
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